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1—1. Abstract 

^■ 

I ' Combining the up-to-date experimental information on i? ^ ttK, ttK* and pK decays, 

^; 

^ . we revisit the decay rates and CP asymmetries of these decays within the framework of 

QCD factorization. Using an infrared finite gluon propagator of Cornwall prescription, 

^ . we find that the time-like annihilation amplitude could contribute a large strong phase, 

Q>^ I while the space-like hard spectator scattering amplitude is real. Numerically, we find that 

tH- ■ all the branching ratios and most of the direct CP violations, except Acp{B -^ K tt^), 

r^ ' agree with the current experimental data with an effective gluon mass rUg ~ 0.5 GeV. 

o ; 

OP ■ Taking the unmatched difference in direct CP violations between B — > tt^K and ir^ K 

O ■ 

^ I decays as a hint of new physics, we perform a model-independent analysis of new physics 

/\^ • contributions with a set of s(l -|- 75)6 Cg) q{l + 75)5 (q=u,d) operators. Detail analyses of 

H : 

5^ , the relative impacts of the operators are presented in five cases. Fitting the twelve decay 

modes, parameter spaces are found generally with nontrivial weak phases. Our results 
may indicate that both strong phase from annihilation amplitude and new weak phase 
from new physics are needed to resolve the ttK puzzle. To further test the new physics 
hypothesis, the mixing-induced CP violations in i? ^ n^Ks and p^Ks are discussed and 
good agreements with the recent experimental data are found. 
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1 Introduction 

With the fruitful running of BABAR and Belle in past decade, plenty of exciting results has 
been produced, which provides a very fertile testing ground for the Standard Model (SM) 
picture of flavor physics and CP violations. Although most of the measurements are in perfect 
agreement with the SM predictions, there still exist some unexplained mismatches. Especially, a 
combination of experimental data on a set of related decays will increase the tension between the 
SM predictions and experimental measurements. At present, there are discrepancies between 
the measurement of several observables in B ^ irK decays and the predications of the SM, the 
so-called "vri^ puzzle" [1] , which have attracted extensive investigations in the SM [21 [3], HI El 
ElEj, as well as with various specific New Physics (NP) scenarios |H]. 

Recently, Belle has measured the direct CP violations B — > Kir decays [9] 

^^ ^ r(50 -^ K--K+) + T{BO -^ K+-K-) ^ ' 

The difference between direct CP violations in charged and neutral modes is 

AA = AcpiB- -^ K^n^) - Acp{B^ -^ A'"7r+) = 0.164 ± 0.037. (3) 

The averages of the current experimental data of BABAR [10], Belle [9], CLEO [H] and 
CDF [12] by the Heavy Flavor Averaging Group (HFAG) [13] are 

Acp{B~ -^ /i"7r°) = 0.050 ± 0.025 , 

Acp(5° -^ A'~7r+) = -0.097 ± 0.012, (4) 

and the difference AA = 0.147 ± 0.028 is established at 5a level. However, within the SM, 
it is generally expected that Acp{B^ -^ n'^K^) and Acp{B^ -^ ir^K") are close to each 
other. For example, the recent theoretical predictions for these two quantities based on the 
QCD factorization approach (QCDF)[T1], the perturbative QCD approach (pQCD)[T5] and the 



soft-coUinear effective theory (SCET) ji(ij read 




f AcpiB- - 
\ Acpim - 


-^ n^K-)QCDF = -3.6% , 

^ QCDF Scenario SA [3] 

-. 7r+K-)QCDF = -4.1% ; 


(5) 


1 AcpiB- - 
[ Acp{B^d - 


pQCD [5 


(6) 


i Acp{B~ - 
1 AcpiBy 


-^ 7rOi^-)sc7ST = (-11 ± 9 ± 11 ± 2)% , 
■^ 7r+K~)scET = (-6 ± 5 ± 6 ± 2)%. 


(7) 



We can see that the present theoretical estimations within the SM are confronted with the 
estabhshed AA. The mismatch may be due to our hmited understanding of the strong dynamics 
in B decays which hinders precise estimations of the SM contributions, but equally possible due 
to new physics effects [TTl [TH] . 

As is known, the annihilation decay of B meson into two light mesons offers interesting 
probes for the dynamical mechanism governing these decays, as well as the exploration of CP 
violation. In most of B meson non-leptonic decays, the annihilation corrections could generate 
some strong phases, which are important for estimating CP violation. However, unlike the 
vertex-type correction amplitude, the calculation of annihilation amplitude always suffers from 
end-point divergence in coUinear factorization approach. In the pQCD approach, such diver- 
gence is regulated by the parton transverse momentum kx at expense of modeling additional 
kr dependence of meson distribution functions [15], and a large strong phase is found. In the 
QCD factorization (QCDF) approach [Tl], to give a conservative estimation, the divergence is 
parameterized by complex parameters, Xa = /q dy/y = ln(mf,/A)(l + p^e*"^-*), with Pa < 1 
and unrestricted (pA, which will sometimes introduce large theoretical uncertainties in the final 
results. In Refs. P, [19], annihilation diagram is studied with SCET and also parameterized 
by a complex amplitude. At present, the dynamical origin of these corrections still remains a 
theoretical challenge. 

In this paper, we will revisit B -^ ttK, ttK* and pK decays within QCDF framework. 
However, we shall quote the infrared finite gluon propagator of Cornwall prescription [20] to 
regulate these divergences in hard-sepctator scattering and annihilation amplitudes. With this 
alternative scheme, we could evaluate both the strength and the strong phase of hard spectator 
and annihilation corrections at the expense of a dynamic gluon mass, which will be fitted in the 



twelve decay modes. It is interesting to note that the infrared finite behavior of gluon propagator 
are not only obtained from solving the well known Schwinger-Dyson equation [SDl EH [22] , but 
also supported by recent Lattice QCD simulations [23]. Numerically, a sizable strength and 
a large strong phase of annihilation corrections are found. Except Acp{B^ -^ K^tt^), our 
predictions for most of the branching ratios and the direct CP asymmetries oi B ^ irK, irK* 
and pK agree with the current experimental data with an effective gluon mass nig = 0.45 ~ 
0.55 GeV. However, we get Acp{B^ -^ K^tt^) = —0.109 ±0.008 which is still in sharp contrast 
to experimental data 0.050 ±0.025. To resolve this mismatch, we perform a model-independent 
analysis of new physics contributions with a set of flavor-changing neutral current (FCNC) 
s(l ± 75)6 ® g(l ± 75)g (q=u,d) operators. To flt the twelve decay modes, parameter spaces 
are found generally with large weak phases. Our results indicate that both strong phase from 
annihilation amplitude and new weak phase from new physics are needed to account for the 
experimental data. 

In Section 2, we revisit B -^ ttK, ttK* and pK decays in the SM with QCDF modifled 
by an infrared flnite gluon propagator for annihilation and spectator scattering kernels. After 
recalculating the hard-spectator scattering and the weak annihilation corrections, we present 
our numerical results and discussions. In Section 3, to flnd resolution to the CP violation 
difference AA, we present analyses of NP operators. Then, using the constrained parameters 
for the operators, we discuss the mixing-induced CP violations in B -^ t^^Ks and p^Ks- 
Section 4 contains our conclusions. Appendix A recapitulates the decay amplitudes for the 
twelve decay modes within the SM [Sj. All the theoretical input parameters are summarized in 
Appendix B. 

2 Revisiting B -^ nK, nK* and pK Decays in the SM 

In the SM, the effective weak Hamiltonian responsible for b ^ s transitions is given as 



r r ^° 



i=3 



C-J^Oj^ ± CsgOsg) 



±h.c., (8) 



where VghV*^ {q = u,c and t) are products of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements [25] , Ci the Wilson coefficients, and Oi the relevant four-quark operators whose explicit 
forms could be found, for example, in Refs. [2112^. 

In recent years, QCDF has been employed extensively to study the B meson non-leptonic 
decays. For example, all of the decay modes considered here have been studied comprehensively 
within the SM in Refs. [21 [3l HI [26]. The relevant decay amplitudes for B — * irK, ttK* and pK 
decays within the QCDF formalism are shown in Appendix A. It is also noted that the frame- 
work contains estimates of some power-suppressed but numerically important contributions, 
such as the annihilation corrections. However, due to the appearance of endpoint divergence, 
these terms usually could not be computed rigorously. In Refs. [21 [3] , to probe their possible 
effects conservatively, the endpoint divergent integrals are treated as signs of infrared sensitive 
contribution and phenomenological parameterized by 



X Ah Jo y 2 

with Pa < 1 and (pA unrestricted. The different scenarios corresponding to different choices of 
Pa and (pA have been thoroughly discussed in Ref . [3] . Although this way of parametrization 
seems reasonable, it is still very worthy to find some alternative schemes to regulate these 
endpoint divergences, as precise as possible, to estimate the strength and the associated strong 
phase in these power suppressed contributions. 

It is interesting to note that recent theoretical and phenomenological studies are now ac- 
cumulating supports for a softer infrared behavior of the gluon propagator [221 [23, [2H]- Fur- 
thermore, an infrared finite dynamical gluon propagator, which is shown to be not divergent as 
fast as ■\, has been successfully applied to the B meson non-leptonic decays [211 [30]. Following 
these studies, in this paper we adopt the gluon propagator derived by Cornwall [20], to regulate 
the endpoint divergent integrals encountered within the QCDF formalism. The infrared finite 
gluon propagator is given by (in Minkowski space) [20] 



where q is the gluon momentum. The corresponding strong coupling constant reads 



Mo Mo 




BiB- 




Figure 1: Feynman diagrams of hard spectator-scattering contributions. 

where /3o = 1 1 — |n/ is the first coefficient of the beta function, and nj the number of active 
flavors. The dynamical gluon mass M'^{q^) is obtained as [20] 



' L ln( -,^ 
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(12) 
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where mg is the effective gluon mass, with a typical value mg = 500 ± 200 MeV, and Aqcd = 
225 MeV. 

2.1 Recalculate the hard-spectator scattering and the annihilation 
contributions 

The next-to-leading order penguin contractions and vertex-type corrections to these decays are 
known free of infrared divergence and well-defined in QCDF [21 El S] , for which we would not 
repeat the calculation and concentrate on the hard-spectator scattering and the annihilation 
contributions. With the infrared finite gluon propagator to deal with the endpoint divergences, 
we will re-calculate the hard spectator and the annihilation corrections in B ^ PP and PV 
decays. The hard spectator scattering Feynman diagrams are shown in Fig. [H where the 
spectator anti-quark goes from the B meson to the final-state Mi meson and the M2 meson is 
emitted from the weak vertex. The longitudinal momentum fraction of the constituent quark 
in the M2(i) meson is denoted by x {y), and C, is the light-cone momentum fraction of the 
light anti-quark in the B meson. To leading power in l/rrib, the hard spectator scattering 
contributions can be expressed as (where x, y ':^ ^ is assumed) 



H,{M,M,) = ^^ [ dxdyd^^^^BiiO^M.ix) 



+ Ml. 



x{y + uj'^{q'^)/^) ^ x{y + uj^{q^)/C) 

(13) 
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Figure 2: Feynman diagrams of weak annihilation contributions. 



for the contributions of operators Q 
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dxdydi "^ > Bi(0'^A/2(^) 



"^A/i {y) Ah 't>rai{y) 



(14) 



for Qi=5,7, and Hi{M^M2) = for g,=6,8- 

In the above Eqs. flTSl) and flT^ . $51 (,^) is the B meson hght-cone distribution amph- 
tude(LCDA), $a/^(x) and (pmiiv) are the twist-2 and the twist-3 LCDAs of hght mesons, re- 
spectively, which are hsted in Appendix B. uj'^{q'^) = Mg{q'^)/M'^, q^ = —Q'^ and Q'^ ~ —^yM"^ 
is the space- hke gluon momentum square in the scattering kernels. The quantities Am^Mz aiid 
BmiM2 collect relevant constants which can be found in Ref. [3]. 

The Feynman diagrams of the weak annihilation topologies are shown in Fig. [21 When both 
Ml and M2 are pseudoscalars, the final decay amplitudes can be expressed as 

X 1 



A\ = Tx I dxdyas{q' 
Jo 

+ ^ 1 .. . . rf rf 0^,(y)0^,(x) 



{xy — Ci;2(g2) + ie){l — xy) {xy — (^^(g^) + ie)x 



xy — uj^{q 



2(n2\ 



le 



$Mi(l/)$M2(a;) 
(15) 
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4 = 0, 
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Jo 
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xy — uj'^{q^) + ie ^ ^ 



(16) 
(17) 
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rf^<P^,{y)^M,{x) 



{xy — uj'^{q'^) + ie){l — xy) ^ 



<^Ha;)0^,(x)<DAf,(2/) 



(18) 



A^ = IT dxdyas{q^)< 



2(1 + a;) 



{xy — uj'^{q'^) + ie)x ^ 



^r-f0™,(y)$Af.(x) 



[xy — uj'^{q'^) + ie)y -^ 



(19) 



where g^ ~ xyM^ is the time-hke gluon momentum square. The "chirally-enhanced" factor 
r*^ is presented in Appendix B. The superscript "z" and "/" refer to the gluon emission from 
initial- and final-state quarks, respectively. The subscript "1", "2", and "3" correspond to three 
possible Dirac structure, with "1" for {V - A) (^ {V - A), "2" for {V -A)®{V + A), and "3" 
for (S — P) ® (5* + P), respectively. When Mi is a vector meson and M2 a pseudoscalar, the 
sign of the second term in A\, the first term in A2, and the second terms in Ag and A^ are 
needed to be changed. When M2 is a vector meson and Mi a pseudoscalar, one only has to 
change the overall sign of Ag. 
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Figure 3: The singularities in integral spaces (left figure) in annihilation contributions and the 
variations of strong coupling constant corresponding to different rrig choices (in unit of GeV). 

As shown by Eqs. (1131) and flT^ of the hard-spectator scattering contributions, the endpoint 
divergences are regulated by the infrared finite form of the gluon propagator. It is easy to 
observe from Eqs. ( TT3l) and ( IT4l) that hard-spectator scattering contributions are real. For 
the annihilation contributions shown by Eqs. (I15p - (TT^ . singularities of the time-like gluon 
propagators at the end-point of integrations (end-point divergence) are moved into integral 
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Table 1: The CP-averaged branching ratios (in units of 10 ^) oi B ^ ttK, ttK* and pK decays 
in SM with different nig (in unit of GeV) are presented in QCDF columns. 



Decay Mode 




QCDF 




Experiment 






nig = 0.3 


nig = 0.7 m 


,g = 0.45 ~ 0.55 


data 


a:- 


-^ n K 


44.4 


16.8 


23.17 ±3.28 


23.1 ±1.0 


b: - 


-^n^K- 


23.4 


9.3 


12.50 ±1.65 


12.9 ±0.6 


X- 


-^ n+K- 


44.7 


16.3 


22.71 ±3.27 


19.4 ±0.6 


Tf.- 




21.2 


7.3 


10.50 ±1.63 


9.9 ±0.6 


BZ- 


-^ TV K 


28.3 


5.2 


8.90 ±1.59 


10.0 ±0.8 


BZ- 


-^ n^K*- 


15.2 


3.4 


5.25 ±0.83 


6.9 ±2.3 


Tf.- 


-^ TT+K*- 


28.7 


5.3 


9.13 ±1.68 


10.6 ±0.9 


Tf.,- 


. n'K*' 


13.4 


1.9 


3.89 ±0.82 


2.4 ±0.7 


B-- 


_T^0 


31.8 


5.6 


10.27 ±1.96 


s.otii 


BZ- 


-^ p^K- 


14.9 


2.5 


4.81 ±0.94 


o O1+0.48 
3.»l_o,46 


52- 


-^ p+K- 


38.6 


8.0 


13.42 ±2.31 


8.61?:? 


iS- 


.p'K' 


21.0 


4.8 


7.53 ±1.25 


c; 4+0-9 

^•^-1.0 



intervals with the infrared finite form of the gluon propagator. Singularities in the integral 
intervals and variations of the effective strong coupling constant are shown in Fig. O It is 
noted that effective strong coupling constant is finite, but rather large in the small q^ region. 
However, there is strong cancellations among the contributions of the small g^ region nearby 
rrtg, which renders the annihilation contribution dominated by g^ > m'^ region associated with 
a large imaginary part. This situation is quite similar to pQCD [IS] where the large imaginary 
part from propagator regulated by kr 

9 p? = Pi 9 tt) -'^'^^{^y^B - ^t)^ (20) 

xyni^Q — kj. + le xym^^ — k^ 

and it is also found the power suppression of these terms relative to the leading contributions 
was not very significant, and important to account for CP violations in i? ^> nK decays. 



Table 2: The direct CP asymmetries ( in unit of 10 ^) oi B ^ -kK, ttK* and pK decays in SM 
with different rUg (in unit of GeV). Other captions are the same as Table 1. 



Decay Mode 




QCDF 




Experiment 






nig = 0.3 


nig = 0.7 m 


,g = 0.45 ~ 0.55 


data 


a:- 


-^ n K 


0.06 


0.19 


0.10 ±0.08 


0.9 ±2.5 


B„-- 


-^TT^K- 


-11.6 


-8.3 


-10.85 ±0.84 


5.0 ±2.5 


X- 


-^ n+K- 


-11.0 


-11.4 


-12.38 ±0.69 


-9.7 ±1.2 


Tf.- 




2.5 


0.1 


1.39 ±0.35 


-14±11 


BZ- 


-^ TV K 


0.3 


-0.0 


0.16 ±0.16 


-11.4 ±6.1 


BZ- 


-^ n^K*- 


-27.0 


-34.1 


-41.20 ±6.69 


4 ±29 


Tf.- 


-^ TT+K*~ 


-27.2 


-47.6 


-47.58 ±8.42 


-10 ±11 


Tf.,- 


. n'K*' 


3.9 


2.1 


4.67 ±1.14 


_q+32 


B-- 


-^ p K 


0.1 


1.2 


0.53 ±0.21 


-12 ±17 


BZ- 


-^ p^K- 


28.1 


49.7 


46.27 ±5.94 


37±11 


52- 


-^ p+K- 


19.3 


31.5 


31.40 ±4.63 


15 ±13 


iS- 


.p'K' 


-4.2 


0.2 


-3.26 ±1.29 


-2 ±29 



2.2 The branching ratios and direct CP asymmetries in the SM 

With the prescriptions for the endpoint divergences, we will present our numerical results of 
branching ratios and CP violations in these decays. Decay amplitudes and input parameters 
are listed in Appendices A and B, respectively. Our results are summarized in Table [T] and 
Table O where the relevant experimental data are also tabled for comparison. 

In Tabled] ([2]), the experimental data column is the up-to-date averages for these branch- 
ing ratios (direct CP violations) by HFAG [13]. It is shown that all the results are in good 
agreements with the experimental data with nig = 0.45 ~ 0.55 GeV. It is also noted that 
the dynamical gluon mass nig = 0.45 ~ 0.55 GeV are also consistent with findings in other 
phenomenal studies of B decays [291 EQ] and the different solutions of SDE [201 EH [22]. The 
phenomenology successes may indicate that the gluon mass, although not a directly measurable 
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quantity, furnishes a regulator for infrared divergences of QCD scattering processes. 
From the CP averaged branching ratios in the fourth column of Table [Tj we get 

'Br{B- -^ n^K- 



Rr 



R„ ^= ~ 
2 



_Br{B- ^n~KO)_ 
Br{B° ^ n+K-y 



1.08 ±0.30, 
1.08 ±0.32, (21) 



BriB^ -^ ttOA'O) 
which agree with the experimental data Re = 1.12 ± 0.10 and /?„ = 0.98 ± 0.09 [T3] . 

Table [2] is our results for direct CP violations. The fourth column is the results estimated 
with rrig = 0.45 ~ 0.55 GeV fixed by branching ratios, where the error-bars are simply due 
to the nig variations. Compared with the experimental data, our results, except Acp{B~ — >■ 
t:^K~), agree with the measurements. For the most significant experimental result among 
the measurements of direct CP violations in the twelve decay modes Acp{B^ — > ti^K') = 
-0.097 ± 0.012 [13], our resuh Acp{B^d -^ tt+R-) = -0.124 ± 0.007 is in good agreement 
with it. As expected in the SM, we find again AcpiB- -^ tt^K') = -0.108 ± 0.008 very 
close to Acp{B^ -^ 7t^K~), which are generally in agreement with the results of Refs. [31E1E] 
listed in Eq. (IH])-(I7]). So, it is very hard to accommodate the measured large difference between 
Acp{B^ -^ Ti^K~) and Acp{B^ — > tx'^K') in the SM with the available approaches for hadron- 
dynamics in B decays. 

Although the problem could be due to hadronic effects unknown so far, the difference 
between Acp{B^ -^ n^K') and Acp{B'^ -^ -k^K~) could be an indication of new sources of 
CP violation beyond the SM [HI EH [32] . 

3 Possible resolution with new (5+P)0(*S'+P) operators 

In this Section we will pursue possible NP solutions model-independently with a set of FCNC 
[S + P)®{S + P) operators. The effects of anomalous tensor and (pseudo-) scalar operators on 
hadronic B decays have attracted many attentions recently [HU [33l [311 [351 [36l ET] . Fo^' example, 
it is shown that they could help to resolve the abnormally large transverse polarizations observed 
in i? ^ (pK* decay, as well as the large Br{B -^ V^*) [36] . 
The general four-quark tensor operators can be expressed as 

O^ = sa^,(l±75)6®ga^"(l±75)g, O^^ = Sia,,,(l ± 75)6^- ® gja'^"(l ± 75)g^ , (22) 

11 



which could be expressed, through the Fierz transformations, as hnear combinations of the 
(pseudo-) scalar operators. In our present case, however, we find that the tensor operators with 



q=u,d give the same contributions to the B^ 



IT 



^K and B^ -^ n^K decays, so that they are 



hardly possible to resolve the direct CP violation difference, because after Fierz transformations, 
O^ and O'^ with q = u,d will give operators like g(l + 75)6 ® s(l + 75)^ which are different 
from s(l + 75)6 ® s(l + 75)5 of the Fierz transforming Of. = S(T^,y(l + 75)6 ® scr^'^(l + 75)5 
for B -^ (f)K* decays. On the other hand, the new operators like s(l + 75)6 ® ^(1 + 75)0' may 
give a possible solution to AA because of their different contributions to the B~ -^ tt^K' and 
B^ -^ 7r+A'~ decays. 

We write the NP effective Hamiltonian for 6 — > s transitions as 



q=u,d 



e^'s 









h.c. 



(23) 



with 0\. and O™ defined by 



Qd 



51 ■= s(l + 75)6®n(l +75)m , 0^8 = Si(l + 75)6j ®Uj(l + 75)ui 



51 - S(l+75)6®t/(l+75)rf, 0^8 = Si(l+75)6j-®t/j(l+75)c?i 



(24) 



where i and j are color indices. The coefficient C^siisS) describes the relative interaction strength 
of the operator O^i/^gw and 5'^ is their possible NP weak phase. Since both the coefficients and 
the weak phase are unknown parameters, for simplicity, we shall only consider their leading 
contributions with the naive factorization(NF) approximation. 



ft(j) 
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5j(i) 



B 



I K. K* 
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IT, p 



(a) 



{h) 



Figure 4: Feynman diagrams contributing to the amplitudes of i? ^ ttK, ttK* and pK decays 
due to the (S + P) (S" + P) operators. 
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The relevant Feynman diagrams of the NP operators are shown in Fig. H] with q = u, d. 
With the NF approximation, it is easy to see that, for the B -^ 7i^K*~ and tt^K*^ decay modes, 
only Fig. H] (a) contributes, for the B -^ n^K^, n~^K~ and pK decay modes, only Fig. H] (b) 
contributes, while both topology structures contribute to the B -^ 7r^K~ and tt^K^ decay 
modes. However, none of them contributes to i? ^ tt'K*^ and 7r^K*~ decays. After some 
simple calculations, these NP contributions to the decay amplitudes of the B — > ttK, ttK* and 
pK decays are obtained as 

Gf 1 



/[NP 
-^B-^ttOR- 



•^B"~^n+R- 



' 71 4 '^*'^** ' ""^^ '^'' ^^ '"^ <^^«) fK , 



V2 4V2 



'tbVts\mB^ 



e^^sgUr^Fr^imDf^ 






BO->nORO 



■ ^P ^ IT/ T/*l 



'tbVts\"''B, 



d'sgj.r^Fi--{ml)fK 






(25) 

(26) 
(27) 

(28) 



where 



/iNP 
•^B-^-K-R*^ 



NP 



•^B-^-k'^R*- 
^NP 

■^BO^nOR*0 



A^^ 



= 0, 

• ^P ^ IT/ T/*l 2 



= 



V2 2V2' *' """ 
Gf 1 



\^6-, g>u _ ^iSl g>d 



■ ^F ^ IT/ T/*l 2 



e^^s g'^^ _ e^^S ^^^ 






r,^<--( 



>^) /^, 



ml) U, 






^NP _ 



/iNP _ 

/iNP _ _ 



9's 



9'i 









9s = C'ss + — Cgi . 



(29) 
(30) 
(31) 
(32) 

(33) 

(34) 

(35) 

(36) 



(37) 
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Comparing the NP amplitudes Eq. (pUj) with Eq. (ITTl) . we expect that these new (pseudo-)scalar 
operators might provide a possible resolution to the direct CP violation difference, which is 
realized in the following numerical analyses. 

3.1 Numerical analyses and discussions of new pseudo-scalar oper- 
ators 

Our analysis consists of five cases with different assumptions for dominance of NP operators, 
namely, 

• Case I: b ^ suu operators O'^^ and O58, 

• Case 11: h -^ sdd operators O'^^ and O^g, 

• Case 111: h —>■ sdd operator O^^ solely, 

• Case IV: only color singlet operators Ogi and Ogi, 

• Case V: all the operators Og^, Ogg, Og^ and Ofg. 

For each case, the corresponding effective Hamiltonian could be read from Eq. (1231) . It could be 
expected that a collection of related decay modes could constrain the relevant NP parameter 
spaces restrictively. 

Our fitting is performed with the experimental data varying randomly within their 2a error- 
bars, while the theoretical uncertainties are obtained by varying the input parameters within 
the regions specified in Appendix B. Our numerical results are summarized in Table [3H5] where 
the assigned uncertainties of our fitting results should be understood at 2a statistical level. 
Illustratively, the constrained NP parameter spaces are shown in Figs. [3HS1 respectively. It is 
noted that, to leading order approximation, both B^ -^ n^K*^ and B'^ -^ -k'^K*' decays do 
not receive these NP contributions, so we perform fitting for the remained ten decay modes. In 
the following, we present numerical analyses subdivided into five cases. 

Case I: 6 ^ suu operators O^^^ and O^g 

We just take into account the contributions of O^^ and O^g in Eq. ( l23l) . i.e. C|^ = C^g = 0. 
In this case, we take the branching ratios of the seven relevant decays B~^7r^K~ , 7r^K*~ , 
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Table 3: The CP-averaged branching ratios (in units of 10~^) in different NP Cases with rUg = 
O.SGeV. The dash means (pseudo-) scalar operators of the Case irrelevant to the corresponding 
decay mode. 

Decay Mode Experiment NP 

data Case I Case II Case III Case IV Case V 



B-^Tx-R 23.1 ±1.0 — 23.0 ±1.0 22.9 ± 0.9 21.5 ±0.3 22.4 ± 0.9 

B-^Tx^K- 12.9 ±0.6 12.1 ±0.4 12.8 ± 0.7 12.7 ±0.6 12.1 ±0.3 12.1 ±0.4 
'^,^Ti+K- 19.4 ±0.6 20.2 ±0.3 — — 20.4 ±0.2 20.1 ±0.4 



X- 


.vrOiT 


9.9 ±0.6 


9.0 ±0.3 


9.9 ±0.6 


10.0 ±0.7 


9.0 ±0.2 


9.1 ±0.4 


B-- 


-^ TrO/i*- 


6.9 ±2.3 


4.2 ±0.2 


4.4 ±0.4 


4.4 ±0.4 


4.3 ±0.3 


4.3 ±0.3 


^.- 




2.4 ±0.7 


3.4 ±0.3 


3.5 ±0.2 


3.5 ±0.2 


3.1 ±0.3 


2.9 ±0.2 



BZ- 


^ p K 


8.0l}i 


— 


B,T- 


-^ pH<~ 


q o 1+0.48 


3.4 ±0.2 


52- 


-^ p+K- 


8.6+?:? 


9.7 ±0.5 


t,- 


.p^K' 


p- 4+0.9 


— 



8.6 ±0.7 8.6 ±0.7 7.4 ± 0.4 7.1 ±0.4 

— — 3.4 ±0.2 3.4 ±0.2 

— — 9.7 ±0.5 9.8 ±0.5 
6.5 ±0.4 6.5 ±0.4 5.5 ± 0.3 5.4 ± 0.4 



p^K~ and B^—*7r^K~, n^K^, tt^K*^, p'^K~ as constraints and leave the direct CP asymme- 
tries as our predictions. The allowed regions of the NP parameters Cg^, C^g and 6g are shown 
in Fig. [5l From which, we find the spaces of Cgi and 6g consist of two parts (dark and gray). 
However, with the gray part, we get Acp{B~ -^ ti^K~) = —0.154 ± 0.038 which conflicts 
with experimental data 0.050 ± 0.025. So, the gray region should be excluded. With the dark 
part of parameter spaces, our prediction Acp{B~ -^ n^K^) = 0.088 ± 0.064 is consistent with 
experimental data. Furthermore, the branching ratios and direct CP asymmetries of the other 
decay modes, listed in the third column of Table and HI agree with experimental data within 
error bars. The constrained parameter space C^^, C^g and 6g are listed in the second column 
of Table El We note that Q^ ^ -Qg ^ -0.04 with S^ ^ 100°, it means the strength of 
color-singlet and color-octet operators are similar, however, such a situation may be hard to be 
generated with a realistic available NP model. 
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Table 4: The direct CP asymmetries ( in unit of 10 ^) oi B ^ irK, ttK* and pK decays. Other 
captions are the same as Table [3] 



Decay Mode 


Experiment 






NP 










data 


Case I 


Case II 


Case III 


Case IV 


Case V 


BZ- 


-^ n K 


0.9 ±2.5 


— 


1.7±2.9 


2.0 ±0.2 


3.9 ±1.0 


3.2 ±1.3 


BZ- 


-^ Tx^K- 


5.0 ±2.5 


8.8 ±6.4 


1.1 ±0.9 


1.2 ±0.9 


2.8 ±5.5 


1.8 ±1.3 


n- 


-^ Tl+K~ 


-9.7 ±1.2 


-5.7 ±4.4 


— 


— 


-10.0 ±0.8 


-9.2 ±1.3 


b:- 


. n^K' 


-14 ±11 


-18.6 ±7.5 


-12.8 ±3.9 


-12.6 ±1.6 


-10.2 ±7.0 


-8.2 ±2.8 


B- - 


-^ TT°K*- 


4 ±29 


4.2 ±19.3 


-8.1 ±3.3 


-8.0 ±3.3 


-4.9 ±19.7 


-13.2 ±4.6 


^^ 




_q+32 


-61.7 ±22.0 


-49.9 ±3.4 


-49.8 ±3.8 


-52.8 ± 24.2 


-47.0 ±6.5 


B- - 


-^ p K 


-12 ±17 


— 


-5.9 ±10.9 


-6.5 ±0.8 


-15.1 ±4.2 


-13.1 ±5.9 


B-- 


-^ p^K- 


37 ±11 


32.8 ±16.5 


— 


— 


48.3 ±3.5 


43.9 ±5.2 


Bd - 


-^ p+K- 


15 ±13 


19.2 ±12.9 


— 


— 


31.9 ±2.7 


28.0 ±4.1 


Bd - 


-^ p^K 


-2 ±29 


— 


-8.1 ±8.1 


-8.5 ±0.9 


-14.9 ±3.0 


-13.5 ±4.4 



Table 5: The numerical results for the parameters C^^, C^^, 5^, Cg^, C^g and (5| in different 
NP Cases. The dashes mean the corresponding operators are neglected in the Case. 

NP para. Case I Case II Case III Case IV Case V 

C^i(xl0-3) -41.6 ±13.4 — — 25.8 ±8.4 -6.7 ±10.5 

C^g(xl0"3) 38.7 ±18.2 — — — 16.0 ±7.1 

(5^ 99.5° ±6.1° — — 107.0° ±11.5° 73.0° ± 23.8° 

C|i(xl0-3) — 23.0 ±5.1 22.8 ±2.3 50.3 ± 12.8 17.5 ±10.1 

C|8(xl0-3) — -0.8 ±13.7 — — 10.5 ±9.4 

(5| — 100.0° ±8.7° 99.3° ±9.2° 106.6° ± 7.3° 114.7° ± 18.6° 



Case II: h -^ sdd operators 0^^ and O^g 

In a large category of NP scenarios with scalar interactions, for example, two-Higgs doublets 
model II, down type fermion Yukawa couplings are enhanced. So, in this case, we evaluate the 
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Figure 5: The allowed regions for the parameters Cg^, C^g and 6^ of Case I 



effects of 0|i and Ofg and neglect Og^ and Og^^. 

As shown by Eqs. (I25])-(I36D, O^^g^ contributes to the decays B^n-R^, tx^K', tx^K*', 
p~K^, TT^K^, TT^K*^, and p^K^. From Table [H one can find that the SM predictions for their 
branching ratios are consistent with the experimental data. So, in this Case, NP weak phase 
Sg would be arbitrary for very small strengths of Cgi and C^g, we thus have to take into 
account both branching ratios and direct CP violations as constraints. The allowed region of 
Cgi, C|g and 6g are shown in Fig. [HI The fitted results are shown in the fourth column of 



Table [HI H] and the third column of Table [5l Interestingly, we note that Cgi = 0.023 ± 0.005, 
Cgg = —0.001 ± 0.013 (consistent with zero) with 6g ~ 100°. It indicates that color-singlet 



operator O^-^ dominates the NP b -^ sdd contributions. Actually, with O*^^ only, we could find 
a solution to the "vr/C puzzle" which will be discussed in next Case. 

Compared with Case I, it is found that |C|J < jC^J ~ IC'ssI- However, we can't conclude 
that 0^^/gN dominates the NP contribution until we consider the two operators simultaneously, 
which will be discussed in coming Case IV and Case V. 

Case III: h -^ sdd operator O'^^ solely 

As the former Case, both branching ratio and direct CP violation are taken as constraints. 
With O^]^ solely, we find a solution to the "vrA' puzzle" with the C|^ and 5'^ allowed region 
shown in Fig. [71 The numerical results are listed in fifth column of Table [3, [Hand fourth column 
of Table [5l respectively. Cg^ and (5^ are found similar to the ones of Case II. It confirms our 
findings in Case II that 0|^ dominates the NP contributions and the contribution of 0|g is 
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Figure 6: The allowed regions for the parameters Cg^, C^g and 6g in Case II with nig = 0.5 GeV. 



negligible. As known, it is easy to generate the situation in many NP scenarios. However, both 
the strength C|^ ~ 0.022 and the new weak phase S'^ ~ 99° normalized to — i-lVJftVJ^I may be 
toughly large for realistic NP models without violating other precise electro-weak measurements. 
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Figure 7: The allowed regions for the parameters Ci, and 6i of Case III 



Case IV: only color-singlet operators Og^ and Og^ 

In order to compare the relative strength of two color singlet operators O^^ and O^^, we 
take them into account at the same time and neglect the other two color-octet ones. Taking the 
branching ratios of the relevant decays as constraints, we find the allowed regions for the NP 
parameters Cg^, 6^, C|^ and 6g, which are shown in Fig. [HI All our predictions for the direct 
CP violations, listed in in sixth column of Table HI agree with experimental data. Especially, 
we note our predictions Acp{B'^n'^K') = 0.028 ± 0.055 and AA = 0.128 ± 0.056 agree with 
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experimental data very well. 

The fifth column of Table [5] is the parameter space obtained for the present Case. We find 
that strength of C|^ in Case IV is larger than the ones in Case II and Case III, because the 
terms of Cg^ and C^^ always have opposite sign in Eqs. fpUl) . fl2Sl) . flHUl) and flH^ . but only 



one of them exists in the other decay modes. It is found that C|^ ^ 2 x C"i ~ 0.05 with 



^51 



6g ~ 6^ ~ 107°, which shows 0^^ dominance. 
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Figure 8: The allowed regions for the parameters Cg^, 5^, C^^ and 5g of Case IV. 



'S8 



Case V: all the operators Og^, O^g, Og^ and Of 

At last, we fit the measured branching ratios and the direct CP violations of all the relevant 
ten decay models with the four operators in Eq. ([21]). Generally the ten CP averaged branching 
ratios are measured with high significants, however, only Acp{B^ -^ ii^K^) is well established 
at 8a level and Acp{B^^7r^K~) by itself is at 2a level only. 

From the fit, the allowed regions for the six NP parameters Cg^, C^g Sg, C|^, C|g and Sg 
shown in Fig. [HI The fitted branching ratios and CP violations are listed in the seventh column 
of Table [3] and HI and the fitted values of the NP parameters are presented in the last column 
of Table O respectively. Since the experimental data are allowed varying randomly within their 
2(7 error-bars, the uncertainties of our fitting results are turned to be quite large. 



We find Qi 



-6.7 ± 10.5) X 10^3 and Qg = (16.0 ± 7.1) x 10"^ with S^ = 73.0° ± 



23.8°, which shift our predication Acp{B^ — > n'^K^) ^ —0.124 in the SM more closer to 
the experimental data —0.097. However, it does not indicate that the b -^ suu operators are 
important for resolving CP violation difference AA, since the sum of their contributions to 
B^ -^ 7i^K~ is quite small due to cancellation among them. For the b — »■ sdd operators, we 
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Figure 9: The allowed regions for the parameters Cg-^^, C^g 5^, C|^, C|g and 6g of Case V. 



get C|i = (17.5 ± 10.1) X 10^3 and C^g = (10.5 ± 9.4) x lO'^ with (5| = 114.7° ± 18.6°. The 
results are consistent with these of Case II and Case III as shown in Table [5], however, due 
to interferences with b -^ suu contributions, the uncertainties are much larger than the two 
former Cases where h -^ suu operators are dropped. Moreover, as shown by Eq. (!26l) . C|g is 
suppressed by 1/Nc in the amplitude of B~ -^ ti^K~ , thus, the dominate status of O'^^ for 
resolving AA is remained. 



3.2 The mixing-induced CP asymmetries in 5 ^ tt'^Ks and B -^ p^Ks 

So far we have discussed the direct CP asymmetries in the these decays with five NP scenarios. 
However, it is naturally to question if we can account for the mixing-induced CP asymmetries in 
B^ -^ TT^Ks and p^Ks decays with these constrained parameter spaces obtained in the former 
subsection. As known, the mixing-induced asymmetries are more suitable for probing new 
physics effects entered via b —^ sqq parton processes than the direct ones, since the former ones 
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Table 6: The mixing-induced CP asymmetries ( in unit of 10 ^) of B^ — > ir^Ks , p^Ks decays. 
Other captions are the same as Table [3] 

Decay Mode Experiment SM NP 

data Case I Case II Case III Case IV Case V 

B^^^n^Ks 38 ±19 77 ±4 45 ±11 56 ± 5 57 ± 3 59 ± 9 62 ± 8 



5° -^ p^Ks 61+^^ 66 ± 3 — 61 ± 6 61 ± 3 56 ± 3 57 ± 4 



could be predicted more accurately in QCDF. Detail discussions for the interesting feature could 

be found in Ref. [3H]- Recently, the measured relative small mixing-induced asymmetry ( with 

large error-bar ) in i?*^ ^ n^Ks has attracted much attention in the literature [38113911^ fTSl I^ . 

The time-dependent CP asymmetries in B^ -^ tt^Ks and p^Ks decays could be written as 

Af{t) = Sf sm{Amdt) - Cj cos(Amdt), (38) 

where —Cf = Acp is the direct CP violation already discussed in former subsection. Sf = A^p 
is the mixing-induced asymmetry 

9Tm A 

A— (5°^/) = -^^ {f = 7r'Ks,p'Ks, Vf = -l) (39) 

^± l^/l 

where A/ = -e-'^'^A^^/A^^ and sin(2/3) = sin(2/3)*Xs = 0.68±0.03[13], since the NP operators 
are irrelevant to B^ — B^ mixing amplitude. 

Using the constrained parameters of the NP operators in Table Eland taking rrig = 0.5 GeV, 
our numerical results are listed in Table for the SM and the five Cases of NP operators. The 
experimental data column is the averages by HFAG [I3]. In the SM, up to doubly Cabibbo 
suppressed amplitudes, one can expect 

Acp ^ 0, Acp = ^ ~ sm{2(3)^Ks = 0-68 ± 0.03 (40) 

for the two decay modes. We get A^^'^iir^Ks) = 0.77 ± 0.04 and A™i?(p°Ks) = 0.66 ± 0.03. 
It is noted that the former is slight larger than sm{2f3)iiJKg which is due to corrections of the 
suppressed amplitudes proportional to VubV*^ as discussed in Ref. [38] 



As shown in Table El 



XU.K 



1 If the old data sin(2/3)vtxs = 0-725 ± 0.037 used, we get AS^ok^, = S^o^^ - sin(2/3)vtKs = 0-05 ± 0.08 and 
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the NP pseudoscalar operators decrease S^o^^ and S'pO^-^ (weaker than former), which seems 
to be favored by the experimental data. 

We note that HFAG has not included the following data yet 



^^]?(5" -^ n^'Ks) = 0.55 ± 0.20 ± 0.03 






0.67 ±0.31 ±0.08 



BABAR 
Belle ^ 



(41) 
(42) 



which are reported very recently at ICHEP08. The average reads A^^p{B^ -^ rc^Ks) = 0.58 ± 
0.17. Again from Table El we find the outputs of all the five Cases with their fitted parameter 
spaces are in good agreements with the new experimental results since the error-bar are still 
large. Taking Case II as example, i.e., assuming NP from b -^ sdd, we present the correlations 
of the direct and the mixing-induced CP asymmetries for B^ -^ t^^Ks and B^ -^ P^Ks decays 
in Fig. [ini where the constrained parameters listed in Tabled] are used. Although all points fall 
in the present experimental error-bars. Fig. fTOl shows interesting correlations between A^*p and 
A^p (S). If the experimental St^oxq shrank to be much lower than sni{2(3)<iKsi the NP Case II 
would give large negative direct CP asymmetry. Similar implication also applies to p^Ks final 
states. 
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Figure 10: Correlation between direct and mixing-induced CP asymmetries for B^ — > n^Ks 
(a) and B^ -^ p^Ks (b) in Case II with nig = 0.5 GeV. The lines are the central value of 
experimental data presented at ICHEP 08. Our plot ranges corresponding the experimental 
error-bars. 



AS'pOjf^ — —0.05 ± 0.07; which agree well with the results in the paper. Considering our different treatments 
of the end-piont divergences, the agreement numerically confirms the observation that the mixing induced CP 
violations are insensitive to strong phases in the decay amplitudes. 
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In summary, assuming NP effects entering B -^ irK, tiK* and pK decays via s{S + P)h ® 
q{S + P)q operators, we have performed fittings for the observables in these decays with a 
model-independent approach. It's found that all the experimental data, especially the direct 
CP violation difference AA, could be accommodated by new b -^ suu or 6 ^ sdd contributions, 
of course by their combination. Assuming the dominance of new h — > suu operators (Case I), we 
find the color-octet operator has the similar strength as the color-singlet one, which is rather 
exotic for electro-weak NP models. However, taking the new h — > sdd operators dominant 
(Cases II and III), we have shown that color-singlet operator s{S + P)h ® d{S + P)d solely can 
provide a resolution to the derivations with a strength about half of 6 ^ suu operators. We 
also have performed fits (Cases IV and V) with both b -^ suu and b -^ sdd contributions to 
infer the their relative size in these decays. It is found that the strength of b ^ sdd is stronger 
than that of 6 ^ suu. In all cases, to account for the experimental deviations from the SM 
predictions for direct CP violations, especially for Acp{B^^7r^K^), new electro-weak phase 
about 100° relative to the SM b — > sqq penguin amplitude is always required. With the fitted 
parameters, we present results for the mixing induced CP asymmetries in 5° -^ t^^Ks and 
p^Ks decays. It is found the NP effects generally reduce 5*^0^^^ and S^oks- However, due to 
the large error-bars, the present experimental data do not further reduce the parameter spaces 
of the NP operators. 

4 Conclusions 

At present, the successful running of the B factories with their detectors BABAR (SLAC) 
and BELLE (KEK) have already taken about 10^ data together at T(45') resonance, and 
have produced plenty of exciting results. Tensions between the experimental data and the 
SM predictions based on different approaches for strong dynamics are accumulated, which 
may be due to our limited understanding of the strong dynamics, but equally possible due 
to NP effects. Motivated by the recent observed AA of the difference in direct CP vio- 
lation between Acp{B^ -^ -k^K^) and Acp{B^ -^ K^ix^) and theoretical issues of end- 
point divergences, strong phases and annihilation contributions in charmless hadronic B de- 
cays, we have revisited the B -^ ttK, nK* and pK decays with an infrared finite form of the 
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gluon propagator supplemented to the QCDF approach. In this way, we can get large strong 
phases from the annihilation contributions, while the hard spectator-scattering amplitudes are 
real. From our numerical analyses, we find that the contributions of the annihilation and 
the hard-spectator topologies are sensitive to the value of the effective gluon mass rrig. With 
rrig = 500 ± 50 MeV, our predictions in the SM agree with the current experimental data well, 
except AcpiB"^ -^ K^tt^). Actually with rrig varying from 300 MeV to 700 MeV, we always 
get AcpiB^ -^ 7r°/i^) ^ Acp{B^ — > K^ir^) as shown in Table [21 which also agree with the 
results in the literature. We conclude that NP effects is required, at least can not be excluded, 
to resolve the discrepancies between the observed AA and the SM expectations. 

With four effective NP b -^ suu and h —>■ sdd operators, we have performed a model- 
independent approach to the discrepancies. Our main conclusions are summarized as: 

• Assuming dominance of 6 — ;► suu operators, the fit gives a quite small center value for 
Acp{B^ -^ K^Tx"^) although consistent with the data within its large error-bar. Moreover, 
the strength of color-octet operator O^g is comparable with color-singlet O^^ which may 
be rather exotic for most NP models. 

• With the h — > sdd operator O^^ solely, the observables m. B ^ nK, irK* and pK decays 
could be well accommodated, since B'^ -^ K^n^ is irrelevant to the b -^ sdd operator 
and it's branching ratio and CP violation agree with the SM prediction very well. 

• Assuming dominance of color-singlet operators Ogi and Og^, it is found that the two 
operators have the similar weak phase with Cg^ ~ \Cg^. 

• For all Cases, to account for the experimental deviations from the SM predictions for 
direct CP violations, especially for AcpiB^^n^K"), new electro- weak phase about 100° 
relative to the SM b —>■ sqq penguin amplitude is always required. 

• With the fitted parameter spaces, the NP operators decrease the mixing-induced CP 
violations in B^ -^ t^^Ks and p^Ks decays, especially that of iv^Ks final states. 

It is reminded that both direct and mixing-induced CP violations have not been well es- 
tablished in most of charmless nonleptonic B decays. Although the difference in direct CP 
asymmetries between Acp{B^ -^ ir^K^) and Acp{B^ -^ A'^vr^) shows some hints of new 
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physics activities, we still need refined measurements of the mixing-induced CP asymmetries 
in the related decays B^ -^ t^^Ks and p'^Ks to confirm or refute the NP hints, since the former 
strongly depends on strong phases in the decay amplitudes while the later not so much and can 
be predicted more precisely. In the coming years, the precision of experimental measurement of 
the observables in these decays will be improved much with LHCb at CERN, which will shrink 
the parameter space and reveal the relative importance of the five Cases studied in this paper. 
Then, the favored Case will deserve detail studies with particular NP models. 
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Appendix A: Decay amplitudes in the SM with QCDF 

The amplitudes for B -^ -kK, ttK* and pK are recapitulated from Ref. [S] 
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p=u,c 



Spu "2 + 2'^3,EW 



(44) 



ylSM 



B0~*7T+K- 



V2^s^ 



Z2 KbV*,A^+K~ 



6pu ai + al + ^4 EW 



SP _ _/9P 



p=u,c 



BO-^nORO 



/ , VpbV*^ < A^ojio 



p=u,c 



-al + -a\ EW - A? + 7^/5; 



3,EW 



C ^ ^ P 

Opu «2 + 2"3,EW 



(45) 



(46) 



where the explicit expressions for the coefficients af = af (M1M2) and Pf = /3f (M1M2) can also 
be found in Ref. [3] . Note that expressions of the hard spectator terms Hi appearing in af and 
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the weak annihilation terms appearing in /3f should be replaced with our recalculated ones. 
The amplitudes oi B ^ nK* and B -^ pK decays could be obtained by setting {ttK) -^ (nK*) 
and (ttK) -^ (pK), respectively. 

Appendix B: Theoretical input parameters 

Bl. Wilson coefficients and CKM matrix elements 

The Wilson coefficients Cj(yu) have been evaluated reliably to next-to- leading logarithmic or- 
der [2HI11J. Their numerical results in the naive dimensional regularization scheme at the scale 



11 = mi, {ph = V^hrrib) are given by 

Ci = 1.074 (1.166), C2 = -0.170 (-0.336), C3 = 0.013 (0.025), 

d = -0.033 (-0.057), C5 = 0.008 (0.011), Cq = -0.038 (-0.076), 

Cj/ae.rr,, = -0.016 (-0.037), Cg/tte.m. = 0.048 (0.095), Cg/ae.m. = -1.204 (-1.321), 

Cio/ae.m. = 0.204 (0.383), Cy^ = -0.297 (-0.360), Cgg = -0.143 (-0.168). (47) 

The values at the scale ph, with nib = 4.80 GeV and Ah = 500 MeV, should be used in the 
calculation of hard-spectator and weak annihilation contributions. 

For the CKM matrix elements, we adopt the Wolfenstein parameterization [13] and choose 
the four parameters A, A, p, and r] as 



A = 0.807 ±0.018, A = 0.2265 ±0.0008, p = 0.141+°;°?9, r/ = 0.343 ± 0.016, (48) 
with p = p (1 — ^) and f] = rj {1 — ^). 

B2. Quark masses and lifetimes 

As for the quark mass, there are two different classes appearing in our calculation. One type is 
the pole quark mass appearing in the evaluation of penguin loop corrections, and denoted by 
rriq. In this paper, we take 

m„ = md = m, = 0, m^ = 1.64 ± 0.09 GeV, mt = 4.80 ± 0.08 GeV. (49) 
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The other one is the current quark mass which appears in the factor ry^ through the equation 
of motion for quarks. This type of quark mass is scale dependent and denoted by niq. Here we 
take [m HE] 



ms{fi)/mg{fi) = 27 A ± 0.4 [M] , m,(2 GeV) = 87 ± 6 MeV [H] , 

mb{mb) = 4.20 ± O.OTGeV [17] , (50) 

where rriq^jj) = (m„ + md)(/i)/2, and the difference between u and d quark is not distinguished. 
As for the hfetimes of B mesons, we take [17] Tb^ = 1.638 ps and Tb^ = 1.530 ps as our 
default input values. 

B3. The decay constants and form factors 

In this paper, we take the decay constants 

fB = (216 ± 22) MeV [50], /b, = (259 ± 32) MeV [50], U = (130.7 ± 0.4) MeV [17], 
/x = (159.8 ± 1.5) MeV [17] /x* = (217 ± 5) MeV [l9], /p = (209 ± 2) MeV [17]. (51) 

and the form factors Il9] 



Fo^^^(O) = 0.258 ± 0.031, F^^^{0) = 0.331 ± 0.041, V^^^' (0) = 0.411 ± 0.033, 
A^~^^* (0) = 0.374 ± 0.034, Af"*^* (0) = 0.292 ± 0.028, ^^^^(0) = 0.323 ± 0.030, 
Aq^p{0) = 0.303 ± 0.029, A^^%0) = 0.242 ± 0.023. (52) 

B4. The LCDAs of mesons and hght-cone projector operators. 

The light-cone projector operators of light pseudoscalar and vector meson in momentum space 
read [HH [3] 



Ml, = ^(^75$p(x)-/iP75^0p(x) 



a/3 



^fy f u^ I \ ^vfv h h 



{M,U = -—[^^vix)- ^^ fc^.fe/ -(^)j,,^ (53) 
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where /ip is defined as mftr^/2, and fp{V) is the decay constant. The chirally-enhanced factor 
appearing in this paper is defined as 



mb(/i)2mg(/i) ' ^ mb(yu)(mg + m^)(/i) 



V/ ^ 2my fy 



^ "^b(/^) fv 

where the quark masses are all running masses defined in the MS scheme which we have given 
in Appendix B2. For the LCDAs of mesons, we use their asymptotic forms [521 [53] 

$p(x) = $y(x) = 6x(l-x), 0p(x) = l, 0^(x) = 3(2x-l). (55) 

As for the B meson wave function, we take the form [51] 

^b{0 = Nb^I - Oexp [ - ( ^/^'' )\^ - iB?] , (56) 

L \Mb — rni)/ J 

where ^b = ^—fntj/Ms, and A''^ is the normalization constant to make sure that Jq d^^siO = 1- 
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